Cavity-Mediated Strong Matter Wave Bistability in a Spin-1 Condensate 
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We study matter wave bistability in a spin-1 Bose-Einstein condensate dispersively coupled to a 
unidirectional ring cavity. A unique feature is that the population exchange among different modes of 
matter fields are accomplished via the spin-exchange collisions. We show that the interplay between 
the atomic spin mixing and the cavity light field can lead to a strong matter wave nonlinearity, 
making matter wave bistability in a cavity at the single-photon level achievable. 
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The macroscopic nonlinear phenomena associated with 
ultracold atoms have become a main stream of research 
interest in the emerging field of atom optics [3]. Mean- 
while they establish intimate connection of the new field 
to other branches of physics, such as nonlinear optics 
and condensed matter physics. One important aspect 
in the field is how to coherently manipulate the nonlin- 
ear behaviors of ultracold atomic ensembles for different 
purposes. 

The ability of an optical cavity to provide feedback be- 
tween input and output light fields can result in the mod- 
ification of the atom-photon interaction in a highly non- 
linear fashion. The exploration of such a nonlinear inter- 
action for applications of both applied and fundamental 
interest has led to many exciting developments, includ- 
ing optical bistability, which was a subject of extensive 
study by the optics community in the 1980's, due mostly 
to the prospect of its use as an optical switch in all- 
optical computers In recent years, rapid technologi- 
cal advancement in cooling, trapping, and condensation 
of neutral atoms has brought new opportunities to cavity 
quantum electrodynamics (QED). A combination of cold 
atoms and large coherence couplings enables single-atom 
trajectories to be monitored in real time with high signal- 
to-noise ratio while allows the vacuum Rabi splitting 
of a single trapped atom to be experimentally observed 

a. 

Instead of a single atom, more recent studies in cavity 
QED focus on cavity systems with a collection of ultra- 
cold atoms d, 0, H, H 3, U I , in which strong coupling 
of ultracold atomic gases to cavity optical field are re- 
alized. This allows us to enter a new regime of cavity 
QED, where a cavity field at the level of a single photon 
can significantly affect the collective motion of the atomic 
samples. This opens up new possibilities to manipulate 
the nonlinear dynamics of ultracold atomic gases with 
cavity-mediated nonlinear interaction. 

So far the works involving in cavity with ultracold 
atomic gases have mainly focused on the interplay be- 
tween the cavity field and the atomic external degrees of 



center-of-mass motion d, i, [13, [lil, [H [H, Ei, [H ■ The 
role of internal spin degrees of ultracold atomic gases in 
the atom-cavity coupling has not yet been seriously ex- 
plored. An intriguing property of spinor Bose-Einstein 
condensate (BEC) is that in addition to the repulsive 
binary collisions, atoms of different spin components 
can couple to each other via spin-exchange interactions, 
which give rise to spin mixing [l6| , a nonlinea r dy namical 
phenomenon under intense theoretical 17 ^ ij, ll9'| and ex- 



perimental investigation 0, m [ial^Tln this work 
we propose a scheme to exploit the atom-cavity coupling 
to control the atomic spinor dynamics in a spinor BEC. 
In contrast to the existing works which focus primarily on 
optical bistability, here we pay particular attention to the 
matter-wave bistability. As we will show, the combina- 
tion of cavity-induced phase shift and the intrinsic spin- 
exchange interaction of a spinor BEC leads to very strong 
matter-wave bistability, providing a new playground for 
exploring the spinor dynamics and cavity nonlinear op- 
tics. 




FIG. 1; Schematic diagram showing the system under considera- 
tion. 



Our model — a spinor BEC with hyperfine spin F ^ 1 
confined in a unidirectional ring cavity — is depicted 
schematically in Fig. [TJ At zero temperature we assume 
single-mode approximation (SMA) that atoms in differ- 
ent spin states can be described by the same spatial wave 
function </!)(r), then each spin component can be associ- 
ated with an annihilation operator Cq (a = ±, 0). A weak 
external magnetic field may be applied to break the de- 
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generacy and provide the quantization axis. The cavity 
is designed in a way that only a single traveling mode 
with frequency tOc, described by an annihilation operator 
a, interacts with the atoms. The cavity is driven by a 
coherent laser field with frequency LOp and amplitude Sp. 
The Hamiltonian under the SMA can be written as 



H = H. 



^c^c+ + cLc_^ — Sc a^a + iSp (a^ — a) 



with Hq describing the dynamics of spinor condensate 
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here Xa is the spin-dependent interaction coefficient [l6[ 
of the condensate. We denote q as the quadratic Zee- 
man shift. Sc = Wp — LOc is the cavity-pump detuning. 
Uo = g^/ (ujp — uja), with g being the dipole coupling 
strength and cua the atomic transition frequency, charac- 
terizes the strength of the atom-photon coupling. We will 
assume that the photon frequency is sufficiently detuned 
away from the atomic transitions so that the atomic up- 
per level can be adiabatically eliminated and the interac- 
tion between photon and atom is essentially of dispersive 
nature. We assume that the photons are tt polarized 
which couple the atoms in the F — 1 ground-state mani- 
fold to the excited manifold with F' = 1. The transition 
selection rule is Amp — 0. However, since the transi- 
tion |F = 1, mp — 0)g \F' ~ 1, m'p — 0)^ is forbidden, 
spin-0 level is not coupled by the photon. For simplic- 
ity, we also assume that the coupling strength between 
the cavity field and spin-± atoms are the same. We will 
treat the leakage of cavity photons phcnomenologically 
by introducing a decay rate k with typical values ~ 1 
MHz. By contrast, the time scale for the atomic spin- 
mixing dynamics is much longer — the measured popu- 
lation oscillation frequency is below 10 Hz for ^^Rb [23j 
and around 50 Hz for ^^Na [24]. This separation of time 
scales allows us to assume that the cavity field always 
follows adiabatically the atomic dynamics 

(1) 



Uo [clc 



The corresponding Heisenberg equations of motion for 
the atomic field operators read 



ic± = c±,Hi 



Uoa)ac±, ICQ — co,Hi 



(2) 



Combining Eqs. ([T]) and in the bad cavity limit one 
can find the effective Hamiltonian H^ff which satisfies 



H, 



Ho 



:^tan-i 



Sc - Uo [clc+ + cLc^ 



(3) 



In the following we adopt a mean-field treatment by 
replacing the operators a and Cq, with the correspond- 
ing C numbers a — (a) and Cq, = \/Na exp [—iOa), 
where JVq and 9^ represent the number and phase of 
the bosonic field for the particles in the spin compo- 
nent a, respectively. We take advantage of the ex- 
istence of two conserved quantities: the total atomic 
number N = + iV_ -|- No and magnetization AI — 
— N^, and simplify our problem into the one de- 
scribed by two variables: the normalized population in 
the spin-0 component x = Nq/N and the relative phase 
6 = 2do — 6*+ —9^. The mean- field counterpart of the 
quantum effective Hamiltonian ([3]) read 



H 



9(1- x) + XaX 

+ U{x), 



1 — X + \/ {1 — xf — Tm? cost 



(4) 



with U {x) = j]^ tan-'^ [Uoil - x) - Sc] /N and m = 
M/N is the atomic polarization. We have defined other 
dimensionless parameters as 

T NXa _ q ^-^ NUo Ep - Sc 

Ki Ki Ki hi h\j 

The equations of motion for x and 9 read 
dx 



dr 
d9_ 
dr 



— 2XaX\/ (1 — x)'^ — -m? sin( 



(5a) 



-2 U 



N 



+ 2Aa 



(1 - x) (1 - 2x) - rn^ 
1 - 2x -I- / ^ ^ — — cos( 

(1 — .t)^ — mn? 



(5b) 



where r = Aci is the dimensionless time. 

From Eq. (j5bp one can see that the cavity modifies the 
atomic dynamics in the same manner as the quadratic 
Zeeman terms g, which will lead to redistribution of 
atomic population among different spin states through 
spin mixing. However, this cavity-induced effective Zee- 
man energy is dependent upon the atomic population 
distribution via Eq. ([T]). It is this inter-dependence of 
the atomic and photonic modes that leads to interesting 
nonlinear dynamics of this coupled system, which will be 
the focus of this work. 

The dynamics of the system can be captured by the 
contour plot of the Hamiltonian iJ, which is intimately 
related to the fixed points {xo,9o) given by the equilib- 
rium solution of Eqs. ([5]). In this work, we only con- 
sider the anti-ferromagnetic atoms (^^Na) with Xa > 0. 
The ferromagnetic case is not qualitatively different. In 
the absence of the cavity field, the equilibrium solutions 
(xq, ^o) have been studied in 
independent solutions of xo 



25|. Besides the phase- 
and xq = 1 — \m\ 
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FIG. 2: (Color online) Phase diagram in the parameter space of 
&c and r\ for different type of solutions: (a) 5 = 0; (b) 6 = it. 
Different regions are differentiated by their colors and are labeled 
with the numbers of corresponding solutions. In the black region, 
no physical phase-dependent solutions can be found. The dimen- 
tionless parameters are estimated to be Aa = 10""^, q ~ 2Aa 
and tJo ~ —5 [2^ . the other parameters are set as m = and 
A'' — 1000. The red dashed line in (a) correspond to rf = 0.8. 



for which the relative phase 6*0 is not well-defined, the 
spinor condensate system supports at most one phase- 
dependent solution with 0o = or tt. The presence of 
the cavity field dramatically changes this property. The 
phase diagram identifying different types of solution is 
mapped out in the parameter space of rj^ and 5c, as shown 
in Fig. [21 We can see that, in certain parameter regime, 
the number of different phase-dependent solutions can be 
more than one, different solution regimes of the coupling 
system can be crossed by varying 5c and . Since these 
two parameters are directly related to the pump laser, 
this means that the dynamical properties of the system 
can be easily manipulated by tuning the intensity or fre- 
quency of the pump laser field. 
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FIG. 3: (Color online) Upper panel: Mean intracavity photon 

1 2 

number \a\ and the normalized spin-0 population X versus cavity- 
pump detuning 5c for the steady-state solutions with rf = 0.8, 
corresponding to the red-dashed line in Fig. (Jfa). The ones that 
represented by the red dotted lines correspond to dynamically un- 
stable solutions. Lower panel: From left to right, the phase-space 
contour plot of H corresponding to different values of 5c marked 
in the upper panel as a, b and c, respectively. 

Here we consider the case with the pump intensity rf' 
fixed, by varying the cavity-pump detuning 5c, the equi- 
librium properties of the system are changed, as shown 
in the red-dashed line in Fig. [HJa). The correspond- 



ing phase-dependent fixed points are derived and the re- 
sults are shown in Fig. [3l The system exhibits typical 
bistable behavior: For certain values of 5c, it supports 
three stationary solutions. A standard linear stability 
analysis shows that in the region with three solutions, 
two of these are dynamically stable and the third one 
is dynamically unstable. Further insights can be gained 
by examing the corresponding contour plot of H (Fig. [31 
lower panel). The unstable fixed points correspond to 
the saddle points in the contour plots. 

The same hysteresis feature can also be identified from 
the mean-field energy diagram by the appearance of a 
swallowtail loop structure (the dots in Fig. 3]). A quan- 
tum calculation involving a direct diagonalization of the 
effective Hamiltonian i?e// confirms a well-known cor- 
respondence between the semiclassical and quantum en- 
ergy levels [27| . namely, in the bistable region the quan- 
tum spectrum (the solid lines in Fig. [5]) exhibits a series 
of anticrossings and when connected, these anticrossings 
form the top segment of the swallowtail structure of the 
corresponding mean-field energy level. (In this example, 
we have, without loss of the essential physics, adopted 
a much smaller system so that the quantum calculation 
can be done within a reasonable computational time.) 




FIG. 4: Quantum and mean-field energy levels with A'^ = 20, 
m = 0.2 and rf — 0.02, the other parameters are the same as 
before. The solid lines are quantum energy levels, the black dots 
refer to the mean-field energy levels with 6 = TV, while the white 
dots refer to those with 6 — 0. 

Let US now return to Fig. [31 As indicated in the up- 
per panel of Fig. [3l both the cavity field and the atoms 
exhibit bistable behavior. The mean cavity photon num- 
ber involved is always less than unity. Remarkably, such 
a small number of photons affect the whole condensate 
and lead to complete population redistribution among 
different internal atomic spin states, which can be read- 
ily observed in experiment. This behavior can be un- 
derstood as following: The collective nature of the con- 
densate greatly enhances the atom-photon coupling such 
that a single photon gives rise to a significant atomic 
phase shift, which in turn strongly modifies the popula- 
tion distribution among the spin states. Bistability re- 
sults from the nonlinear feedback between photons and 
atoms. 

It is interesting to compare our study with the experi- 
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mental work of Refs. 11|, where the motional degrees 
of freedom of ultracold atomic gases represent the source 
of nolinearity affecting light-atom interactions. In their 
system the atomic zero-momentum mode and two side 
modes with momentum ±2hk (fc is the cavity light wave 
vector) dominate the dynamics. The nonlinearity orig- 
inate from the coherence between these modes which is 
induced by the coupling provided by the standing-wave 
cavity field [2§|. Such a coupling is neither strong as 
it is provided by the weak cavity field, nor resonant as 
there is a detuning of Ahujrec with tOrec = hk^ /2m the 
atomic recoil frequency. In these systems, optical bista- 
bility of low photon numbers is possible and is indeed 
observed. In principle, atomic population in different 
motional states should also exhibit bistability. However, 
such matter-wave bistability will be difficult to observe 
as it is not easy to measure atomic population of different 
momentum states inside a cavity in real time. Further- 
more, due to the inefficient coupling between atomic mo- 
mentum modes as we just mentioned, the matter-wave 
bistability is very weak since most of the atomic pop- 
ulation will remain in the zero-momentum state. In a 
recent work [ISJ where this system is theoretically exam- 
ined, it is found that bistable behavior involving tens of 
photons can only transfer about 20% of the total atomic 
population out of the zero-momentum state (see Fig. 1 



of Ref. [15|). 



In contrast, in the model we considered here, the 
coherence between the internal atomic spin states af- 
fecting atom-light interaction is induced by the intrin- 
sic spin-exchange interaction which represents a matter- 
wave analog of the four-wave mixing in nonlinear optics. 
An immediate advantage is that it can be independently 
tuned with respect to the cavity field, thereby dramat- 
ically increasing the chance of large population change 
among spin states at a low cavity field. Consequently, 
our system can exhibit very strong matter-wave bistabil- 
ity. This is indeed confirmed by our detailed calculations. 

In summary, we have studied the mutual interaction of 
a spinor condensate with a single-mode cavity field. We 
show that the coupled cavity-spinor condensate system 
can display simultaneously strong optical bistability at 
the single-photon level and strong matter-wave bistabil- 
ity involving a whole condensate with macroscopic num- 
ber of atoms. This opens up new opportunities to explore 
a diversity of new phenomena in cavity nonlinear optics 
with low photon numbers and many-body physics with 
quantum gases. Before ending, we note that the conden- 
sate depletion may become significant in the long time 
scale as the quantum fluctuations of the cavity can in- 
troduce excess noise to the condensate system [2^. A 
more careful treatment taking these effects into proper 
account will be left for further investigation. 
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